Cross-linker unbinding and self-similarity in bundled cytoskeletal networks 
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The macromechanical properties of purely bundled in vitro actin networks are not only determined 
by the micromechanical properties of individual bundles but also by molecular unbinding events 
of the actin binding protein (ABP) fascin. Under high mechanical load the network elasticity 
depends on the forced unbinding of individual ABPs in a rate dependent manner. Cross-linker 
unbinding in combination with the structural self-similarity of the network enables the introduction 
of a concentration/time superposition principle - broadening the mechanically accessible frequency 
range over 8 orders of magnitude. 

PACS numbers: 87.15.-v, 87.15. La, 83.60.Df 



The cytoskeleton is a highly complex, dense and het- 
erogeneous network which offers the cell the crucially 
needed dynamic adaptability of its mechanical behav- 
ior [l[ . For this purpose each cell exploits the whole range 
of accessible building blocks such as actin filaments, mi- 
crotubules and intermediate filaments as well as their as- 
sociated binding proteins to form localized cytoskeletal 
structures adapted to their special needs 0]. An impor- 
tant feature the cytoskeleton has to fulfill is the ability 
to withstand deformations on short time scales while still 
allowing the cell to adjust and rearrange for processes 
taking place on longer time scales. 

The semi-flexible nature of actin filaments themselves 
results in a separation of length scales and therefore time 
scales as already observed for in vitro entangled actin 
solutions [H, H[ . The addition of actin-binding proteins 
(ABPs) alters the structure of entangled solutions de- 
pending on the molecular structure of the cross-linking 
molecule used into isotropically cross-linked networks [j| , 
purely bundled networks @, 0] or heterogenous compos- 
ite phases [8, 9] . Concomitant with the structural change 
new length and time scales are introduced. While the 
microstructure of such complex reconstituted actin net- 
works mainly determines the elastic response 13, 11, 
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the length dependence of the persistence length of indi- 
vidual bundles has to be considered as well M [3, 
Additionally, the typical lifetime of an actin/cross-linkcr 
bond seems to influence the frequency behavior of cross- 
linked actin networks. 

Here we show, that the viscoelastic behavior of bun- 
dled actin/fascin networks crucially depends on molecu- 
lar unbinding events of actin/fascin bonds. On interme- 
diate time scales comparable to the cross-linker off-rate 
the non-linear mechanical response is highly sensitive to 
the applied strain rate. On the same time scale, the 
frequency dependence of the viscoelastic moduli shows 
a generic behavior which is a signature of unbinding 
events and allows creating a master curve. Such a cross- 
linker concentration/time superposition enables us to de- 
termine the viscoelastic properties over almost 8 decades 



in rescaled frequency. 

G-actin is obtained from rabbit skeletal muscle and 
stored in lyophilizcd form at 



15 . For measure- 



-21 °C 

ments the lyophilized actin is dissolved in deionized water 
and dialyzed against G-Buffer (2 mM Tris, 0.2 mM ATP, 
0.2 mM CaCl 2 , 0.2 mM DTT and 0.005 % NaN 3 , pH 8) 
at 4 °C. The G-actin solutions are kept at 4 °C and used 
within seven days of preparation. The average length 
of the actin filaments is controlled to 21 um using gel - 
solin obtained from bovine plasma serum following [lq . 
Recombinant human fascin (55 kD) was prepared by a 
modification of the method of [13] as described by [l8j |. 
In the experiments the molar ratio R between fascin and 
actin, R = Cf/c a , is varied over almost three decades. 

The viscoelastic response of actin/fascin- networks 
is determined by measuring the frequency-dependent 
viscoelastic moduli G'(u>) and G"(ui) with a stress- 
controlled rheometer (Physica MCR 301, Anton Paar, 
Graz, Austria) within a frequency range of three decades. 
Approximately 520 ul sample volume are loaded within 
1 min into the rheometer using a 50 mm plate-plate ge- 
ometry with 160 pm plate separation. To ensure linear 
response small torques (« 0.5 pNm) are applied. Actin 
polymerization is carried out in situ, measurements are 
taken 60 min after the polymerization was initiated. 

Above a critical concentration of the actin binding 
protein (ABP) fascin, actin filaments are organized into 
a purely bundled network. The degree of cross-linking 
among individual bundles depends on the fascin concen- 
tration Q. The non- linear elasticity of fascin networks 
in the purely bundled phase is probed best by applying 
a constant shear rate dj/dt. The resulting stress a is 
reported and from the smoothed (7(7) relation the differ- 
ential modulus K = da/dj is calculated. With increas- 
ing cross-linker concentration R the non-linear response 
exhibits a transition from a strain-hardening regime to 
strain- weakening as depicted in Fig. |TJl. As with increas- 
ing R the bundle stiffness as well as the connectivity of 
individual bundles is enhanced, at high fascin concen- 
trations the macroscopic deformation will lead to forced 
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FIG. 1: Differential modulus K normalized by its value in 
the linear regime as a function of strain. The non-linear 
response can be tuned from hardening to weakening at a 
fixed strain rate d*y/dt = 6.25%s~ 1 by varying the fascin 
concentration (R — 0.02,0.05,0.1,0.2,0.5) as shown in a) as 
well as for fixed R = 0.1 tuning the strain rate (dy/dt — 
0.125, 1.25, 6.25, 12.5%s _1 ) as depicted in b). 
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FIG. 2: Differential modulus K normalized by its value in 
the linear regime as a function of strain. A weakly bun- 
dled network (R = 0.05) is examined at different strain 
rates (d-y/dt = 0.125, 0.625, 1.25, 6.25, 12.5 %s _1 ). The in- 
set shows the maximum stress a max vs. the strain rate. 



unbinding of interconnecting fascin molecules even before 
non-linear behavior can be evoked. 

On the other hand the same effect can be achieved by 
applying different strain rates to a bundle network with 
fixed microstructurc {R = 0.1) as shown in Fig. lb. 
Varying the strain rate by two orders of magnitude leads 
to a continuous change in the degree of strain-hardening. 
This underlines the non-universal behavior of the non- 
linear response of semi-flexible polymer networks as also 
reported for purely entangled actin solutions [HI and 
isotropically cross-linked networks When the bundle 
network is sheared with very low strain rates, detaching 
of fascin molecules is very likely on the timescale of the 
experiment (w 100 s in the case of d^/dt = 0.125 %.s _1 ). 
Therefore, the remaining number of connection points 
between fascin bundles might become too low to evoke 
strain hardening as in the case of higher shear rates. This 
result suggests that the elastic response of purely bundled 
actin networks might be equally determined by the cross- 
linker concentration and the chosen time scale of the ex- 
periment. The non-static nature of actin/fascin bonds 
introduces the cross-linker off-rate as an important time 
scale. For a single molecular bond a linear increase of the 
maximal loading force with the logarithm of the loading 
rate, f ma x ~ ln(df/dt), would be expected 2(| 21 1 . 

To this end we conducted shear rate experiments 
in such a way, that the non-linear reponse was com- 
pletely tuned from strain-hardening to weakening. To 
further investigate the dependence of fmax on the load- 
ing rate, a weakly bundled network is chosen where 
the density of bundle interconnection points is low 
(R = 0.05). For this network type strain- hardening 
is observed for all strain rates applied (FIG. [2|). The 
maximum strain the bundle network can endure with- 
out being irreversibly damaged is independent from the 
strain rate, j m ax ~ (22 ± 1) %. The maximum stress 



&max = c(7mai) the bundle network can withstand de- 
pends loarithmically on the strain rate &y /dt as depicted 
in the inset of FIG.O The maximum force a single fascin 
molecule connecting two bundles can hold is given by 
fmax ~ Cmaxlc whereas l c denotes the average dis- 
tance of two neighboring bundle interconnection points 
beeing fixed for a given R. For the data presented here 
the proportionality df/dt ~ dj/dt holds. Therefore the 
Bell-prediction is reproduced surprisingly well - indicat- 
ing that forced unbinding of actin/fascin bonds gives rise 
to the observed strain rate dependence. As the macro- 
scopic deformation is transmitted in a non-affine way to 
the bundle crosslinking points jf| [22|, not all intercon- 
necting fascin molecules will be loaded with the same 
force during a strain rate experiment. Surprisingly, the 
Bell prediction is still fulfilled. 

Unbinding of fascin molecules between single bundles 
should also be observable in the low frequency behavior 
of the viscoelastic moduli. Indeed, a pronounced mini- 
mum in G" \lo) can be found in the linear response of such 
networks. From the frequency spectra we roughly esti- 
mate the time scale for these events to be on the order of 
several seconds which in general matches the typical off- 
rates measured for various actin-binding proteins [23l. I24I] . 
The measured frequency spectra of actin/fascin networks 
show two subpopulations in dependence on R: Inter- 
connected bundle networks exhibit a plateau-like storage 
modulus G'(u>) while the loss modulus G"(uj) reveals a 
well-defined minimum which shifts to higher frequencies 
with increasing fascin concentration. Networks below the 
bundling transition (low R) 7 however, feature frequency 
spectra that do not exhibit a minimum in G"(uj) but 
nevertheless highly resemble each other. 

To obtain appropriate parameters for a putative gener- 
alization process an impartial criterion is needed. There- 
fore for each frequency spectrum the minimum position 
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in the loss modulus (uj*\G"*(u))) together with the cor- 
responding value of the storage modulus is determined. 
Normalizing the frequency spectra by these values re- 
sults in a master curve which shows two different regimes 
corresponding to distinct network types. Therefore, the 
master curve is discussed best in two parts: Firstly, at 
low rescaled frequencies the master curve corresponds to 
networks in the purely bundled phase (FIG.EK, rescaled 
G'(lu) shifted up for clarity). This regime contains a pro- 
nounced plateau region in the storage modulus G'(uj) 
accompanied by a clear minimum in the loss modulus 
G"(w). Secondly, networks before the bundling thresh- 
old also show a common frequency dependence similar to 
a purely entangled actin solution (FIG.EJd). 

For an actin concentration of c a = 9.5 pM the thresh- 
old ratio R* was determined to be R* ~ 0.01 @. Indeed, 
the scaling parameters obtained for the master curve con- 
struction follow a conjoint power law G'* ~ (u)*) x , 
G"* ~ (uj*) x in the purely bundled phase, while drop- 
ping off this relation for values R < R* (see inset of 
FIG. [3b). For the pure bundle phase a linear relation 
between the scaling parameters, x = 1, is obtained. 
This is an indication of self-similarity as also observed 
in other soft matter systems (25[. Fascin networks with 
R < R* can also be rescaled with a common power 
law between the scaling parameters, however x ~ 1/2. 
As in this regime the tube model can be applied to de- 
scribe the mechanical response 0], the addition of few 
fascin molecules seems to only slightly modify the net- 
work properties which does not give a self-similar struc- 
ture. 

The corresponding part of the master curve, 
(jj/lj* > 10 3 , shows a regime ~ a> ' 5 in both viscoelastic 
moduli right before the crossing-over, Lo cross ~ l/r e , 
where r e denotes the entanglement time. While the 
molecular origin of this behavior is still unclear, it is 
expected for cross-linked scmi-flcxiblc polymer networks 
in an intermediate frequency regime resulting from ten- 
sions in the network [4[. It is also reported for entan- 
gled actin networks from two-point micro-rheological ex- 
periments [26[ and interpreted to be due to the diffu- 
sive dissipation of long-wavelength longitudinal fluctu- 
ations. Only at very short time scales t << r e , and 
therefore very high frequencies beyond the crossing-over 
of G'(ui) and G"(u>), i.e. around 10 kHz, a scaling 
~ oj - 75 due to undisturbed relaxations along single fil- 
aments would be expected 0, [13, HH . Rcscaling the loss 
factor tan(<5) = G" (ui) / G' (u) (supplementary informa- 
tion) shows unambiguously that the w 5 regime observed 
for actin/fascin networks is indeed an analytical power 
law as the Kramers-Kronig relations are fulfilled. Possi- 
bly, the frequency behavior at intermediate times might 
be strongly dependent on the network architecture as a 
regime ~ cj 0,75 has been described for actin networks 
in a composite phase, as obtained by the cross-linker 
scruin 1291. 
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FIG. 3: Master curve for actin/fascin networks a) in the bun- 
dle phase (R = 0.5,0.2,0.1,0.05,0.02,0.01) and b) before 
the bundling transition (R = 0.01, 0.005, 0.002, 0.001, 0) as 
described in the text. Closed symbols denote G' , open sym- 
bols denote G" . The inset shows the dependence of G'* and 
G"* on u* (full line: ~ (a;*) 1 , dotted line: ~ (oj*) 1/2 ). 



Thus the mesoscopic network structure does not only 
determine the plateau modulus but also the frequency 
behavior of semi-flcxiblc polymer networks. Despite the 
fact that for composite networks a master curve descrip- 
tion was reported [29^ , this is not a generic feature for 
cross-linked actin networks. Purely isotropically cross- 
linked networks [j| do not allow a master curve construc- 
tion based on a concentration/time superposition. There, 
the network elasticity is enhanced by simply decreasing 
the cross-linker distance l c - without any overall change 
in the network structure. Thus self-similarity does not 
apply for isotropically cross-linked networks. The mas- 
ter curve construction applied here is only possible if one 
intrinsic parameter determines structure and mechanical 
behavior at the same time, which is the case for purely 
bundled actin/fascin networks. Raising the fascin con- 
centration R "magnifies" all system properties includ- 
ing the mechanical parameters of its constituting bun- 
dles @, [IH ■ This creates a coarsening, self-similar net- 
work that can be described by the principle of cross-linker 
concentration/time superposition introduced here. 
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Still, thermal and enforced unbinding of interconnect- 
ing fascin molecules could lead to subtle changes in the 
network structure. The increase of G"(cu) at frequen- 
cies lower than to* is quite weak, which stands in marked 
contrast to rigor HMM-networks where a steep uprise of 
G"(u) below the minimum position is reported. This 
difference is attributable to the fact that the mechanical 
response of an isotropically cross-linked network is dom- 
inated by the number of HMM-molecules bound, corre- 
sponding to l c . Thus unbinding results in a structural 
and mechanical transition from a cross-linked network to 
an entangled solution. In contrast, for a purely bundled 
network the linear response is dominated by the bundle 
properties and the non-afhnc deformation mode. Only at 
large time scales t >> 1/lo* it will be slightly affected by 
the cross-linking degree. 

In summary, the binding kinetics of the cross-linking 
molecules and the network structure determine the net- 
work elasticity at high deformations. Forced unbind- 
ing of cross-links under mechanical load can tune the 
non-affine response on intermediate times. Moreover, 
the coarsening, self-similar structure of purely bundled 
actin/fascin networks allows a generalization of their fre- 
quency behavior. The concentration/time superposition 
introduced here broadens the measurable frequency win- 
dow over 8 orders of magnitude simply by altering the 
typical network length scales tuning R - fully equivalent 
to the principle of temperature/time superposition used 
for classical flexible polymer networks. Thus, the chosen 
time scale crucially affects the mechanical network re- 
sponse as distinct length scales dominate the viscoelastic 
behavior. As a consequence a generic broad distribution 
of relaxation times has to be considered describing the 
frequency behavior of cross-linked semi-flexible polymer 
networks and living cells 3^, 3l| . The self-similar in vitro 
system presented here provides a benchmark to address 
the challenging question of how single filament and bun- 
dle properties, network structure and binding kinetics of 
ABPs interplay to provide the high adaptability, which 
is known for the cytoskeleton of living cells. Therefore, 
the described effects may be more generic for an under- 
standing of adaptable biomaterials. 
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Supplementary information 
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FIG. 4: Generalized loss factor curve as described in 
the text. The unshifted curves are depicted in the inset 
(R = 0,0.001,0.002,0.005,0.01,0.02,0.05,0.1,0.2,0.5). 

From a single macrorhcological experiment it is dif- 
ficult to extract the crossing-over time T e with a decent 
accuracy. The master curve presented in this work would 



in principle offer enough reliable data points - albeit nor- 
malizing both moduli to the same absolute value cancels 
out the desired information. In fact, shifting up again 
the rescaled storage modulus might give the wrong im- 
pression that both moduli do not cross over at all in the 
frequency range probed. Therefore, it is helpful to calcu- 
late the loss factor tan (S) = G" (ui) / G' (ui) . The shape of 
the loss factor curves drastically changes at R* — 0.01 
(see inset of supplementary FIG. 2]). This again is in 
agreement with the bundling transition reported before. 

Moreover, these curves can also be generalized (FIG. 01 
requiring only one single rescaling parameter, namely the 
same to* as used for the master curves shown in FIG. [3] 
For a rescaled frequency range of 10 3 — 10 4 the loss 
factor master curve exhibits a plateau. In this regime 
the viscoelastic moduli show the same frequency depen- 
dence and have comparable absolute values. Thus the 
Kramers-Kronig relations are fulfilled implying an ana- 
lytical power law behavior in an intermediate asymptotic 
regime. For even higher rescaled frequencies the loss fac- 
tor finally reaches values larger than 1 , indicating that a 
crossing-over to a viscous dominated regime, a; > l/r e , 
does exist but is not sufficiently accessible with macrorhc- 
ological methods - even not by the application of the 
cross-linker/time superposition. 



